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Abstract. The light propagation through a stack of two media with different refractive indices, which are 
aligned according to different quasiperiodic sequences determined by metallic means, is studied using the 
transfer matrix method. The focus lies on the investigation of the influence of the underlying quasiperiodic 
sequence as well as the dependence of the transmission on the frequency, the incidence angle of the light 
wave and different ratios of the refractive indices. In contrast to a periodically aligned stack we find 
complete transmission for the quasiperiodic systems for a wide range of different refraction indices for 
small incidence angles. Additional bands of moderate transmission occur for frequencies in the range of 
the photonic band gaps of the periodic system. Further, for fixed indices of refraction we find a range of 
almost perfect transmission for angles close to the angle of total reflection, which is caused by the bending 
of photonic transmission bands towards higher frequencies for increasing incidence angles. Comparing with 
the results of a periodic stack the quasiperiodicity seems to have only an influence in the region around 
the midgap frequency of a periodic stack. 



1 Introduction 

Since the discovery of quasicrystals by Shechtman et al. in 
1984 [1 , many efforts have been conducted to understand 
the physical properties of these aperiodic materials, which 
possess a long range order without having a translational 
symmetry. Hence, quasicrystals are regarded to have a 
degree of order intermediate between crystals and disor- 
dered systems. Quasicrystalline systems have been exten- 
sively studied e.g. with respect to their structure, which 
shows uncommon rotational symmetries |1|2|3] , and their 
electronic states, which were found to show a Cantor-set 
spectrum in one dimension |4)5|6|7j ; but also phonons and 
magnetic properties of these materials have been investi- 
gated |8I9I10I11I12| . 

Further, the photonic properties are of special inter- 
est because the complex symmetries in quasicrystals make 
them suitable for the application in several optical de- 
vices such as single-mode light-emitting diodes, polariza- 
tion switching and microelectronic devices that are based 
on photons rather than on electrons, which potentially 
can be the electromagnetic analogue to semiconductors 
113114115116117] . 

For the theoretical study of quasiperiodic systems one 
often applies the concept of aperiodic mathematical se- 
quences / tilings. Especially, the photonic properties of 
one-dimensional systems have been extensively analyzed 
with this approach, not only because transfer matrix meth- 
ods can be applied but also because such systems can be 
relatively easily produced in reality. There are examples 
based upon the Fibonacci sequence, the Thue-Morse se- 
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quence or Cantor sequences |14|18|19(20] , and also systems 
for negative refractive indices have been studied [21 22 23 . 
Further, one can find a close resemblance of the theoretical 
and the experimental results |24|25j . 

However, most of the investigations do not deal with 
the angle of incidence, but rather assume that the inci- 
dent light wave is perpendicular to the layers of the stack, 
and do not investigate the influence of the structure of 
the quasiperiodic sequence on the transmission. In this 
paper we discuss the photonic properties of quasiperiodic 
systems consisting of a sequence of layers made of two dif- 
ferent refractive indices ha and n# (cp. Figure [TJ, where 
the configuration of the layers is given according to the 
so-called metallic-mean sequences. In particular, we focus 
on the transmission of the light wave in dependence on 
the incidence angle $ and investigate the influence of the 
underlying construction rule for different metallic-mean 
quasicrystals on the photonic properties as well. We find 
that additional transmission bands occur for frequencies 
corresponding to the photonic band gap of the periodic 
system. Further, for small incident angles we obtain al- 
most complete transmission for the quasiperiodic systems 
for a wide range of ratios of refractive indices ua/tlb in 
contrast to a system with periodically stacked layers, and 
for angles close to the angle of total reflection there is also 
a region of relatively high transmission. 

The configuration of the layers of these so-called metal- 
lic-mean quasicrystals can be constructed by the inflation 
rule 



A ABA a ~ x 
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Fig. 1. Transmission of light through a quasiperiodic 
stack consisting of the materials A and B with differ- 
ent indices of refraction. The light wave is incident 
on the surface with an angle Parts of the light wave are 
transmitted through the stack (E^) and other parts are 
reflected (E?). 

iterated ra times starting from the symbol B. For instance, 
this yields for the parameters a = 1,2,3 the approximants 
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We refer to the resulting sequence after m iterations as 
the rath order approximant with the length f m given 
by the recursive equation f rn = f rn -2 + a/m-i and fo = 
fi = 1. Depending on the parameter a the inflation rule 
generates different metallic means, i.e., the lengths of two 
successive sequences satisfy the relation 



frr 



lim 



(3) 



where r(a) is an irrational number with the continued 
fraction representation r(a) = [a] = [a, a, a, ...]. Also, the 
number of layers A and number of layers B in 
an approximant are related by these metallic means 
according to 

lim ||=T(a). (4) 

For example, a = 1 yields the well-known Fibonacci se- 
quence related to the golden mean ta u = [1] = (1 + a/5)/2, 
the case a = 2 corresponds to the octonacci sequence with 
the silver mean t Ag = [2] = 1 + >/2, and for a = 3 one 
obtains the bronze mean tb z = [3] = (3 + y/l3)/2 |2Q|26j . 
In general the relation r(a) = (a + \J a 2 + 4)/2 holds [27] . 

Due to the recursive inflation rule , these quasiperi- 
odic chains possess a hierarchical structure, which is more 
clearly visible by using the alternative construction rule 



na na na (na \a — 1 

u ra — u m-l u m-2 V u m-1/ 



(5) 



yielding the same quasiperiodic sequences for ra > 2 set- 
ting Cq = B and C± = A. Likewise, there is a correspond- 
ing expression for the recursive construction of the transfer 
matrix as given in Sec. [2] 



The outline of this paper is as follows: In Sec. [2] we 
give an introduction to the transfer matrix method used 
for the calculations and in Sec. [3] we present our results 
and discuss them. This is followed by a brief conclusion in 
Seed 



2 Transfer Matrix Method 

The propagation of light through a layered system as shown 
in Figure [l] is commonly investigated by means of transfer 
matrix methods [18 1 9|2l(23] , Regarding the geometry we 
consider the propagation of linearly polarized light with 
the electric field perpendicular to the plane of the light 
path (transverse electric waves) [28] . 

Within one layer the electric field E 
described as a superposition of a right and a left traveling 
plane wave 



Ee y can be 



E 



■E{e 



ikjX — iojt 



(6) 



where kj = njk denotes the wave number of the light in 
the jth layer made either of medium A or B with refrac- 
tive index ra,, k denotes the wave number in the vacuum 
and uo denotes the frequency of the light. Note, that each 
amplitude E 3 V and E\ actually represents the resultant of 
all possible waves traveling in the particular direction. 

The propagation of the waves through the stack is 
given on the one hand by their refraction at the inter- 
faces between the layers and on the other hand by their 
propagation through the layers. Addressing the first point, 
the boundary conditions at the interfaces between the lay- 
ers require the tangential component of the electric field 
E and the magnetic field H = x E to be continuous 
across the boundaries. This yields for an incidence angle 
?9 a , the emergence angle and the corresponding indices 
of refraction n a and np the relations (cp. Figure [2| 



El + E{ 
(E 3 X — E J )n a cos $ a 



E^ 1 - 
(Et 1 



El +l 

7.3+1 



(?) 



£;/ +i )n^cos^ . (8) 



Applying the notation used by Kohmoto et al. [18 19 with 
the variables 

E+ = E r + Ei and E_ = -i (E r - E\) , (9) 

the transfer of the light wave from medium a to medium 
j3 can be described by 



E_ 



- T a ( E+ 



j + 1 



with the interface matrix 



T n R — To 



l a(3 — -L(3 a ~ I q n p cos ^ 
v n a cos 'dry 



l o 



(10) 



(11) 



Thereby, the incidence and emergence angles are related 
by Snell's law [18] 



Sin^ a Uf3 

sin^ n a 



(12) 
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Fig. 2. Electric and magnetic field vectors E and H at the 
interface between layer j and j + 1 with different indices 
of refraction n a and np. 



On the other hand, the propagation of the light waves 
within the layers results in a phase difference, which is 
comprised by the propagation matrix 



cos ip 1 sin (p 7 
- sin cp 7 cos ip 7 



(13) 



The phase difference for the wave length A = 2ii/k of the 
light and a layer thickness <i 7 amounts to (cp. Figure [3| 



<p 7 = k[n 7 (0P + PQ) - n a OR] 



47T 



<i 7 n 7 cos$ 7 . (14) 

A 

Note, that the paper by Kohmoto et al. [18 and some 
other related papers as e.g. [20 25 only state the phase 
shift, which the waves undergo by traversing the stack, 
and not the phase difference between adjacent waves. 
The transfer matrix 



(15) 



N 



of the system is then given as a combination of the dif- 
ferent interface and propagation matrices according to a 
certain quasiperiodic sequence. Thereby, the amplitude of 
the incident light is denoted by E®, of the reflected light 
by E® and of the transmitted light by E^ (cp. Figure [l]) 

For instance for the Fibonacci sequence ABAAB . 
one obtains 



M = TcaTaTabTbTbaTaTaTabTb 
• • • T b T B aTac • 

In general the recursive equation 

mZ = aC_ 1 aC 2 {aCi}°" 1 



(16) 



(17) 



is applicable for m > 2 with Mq = TabPbPba and 
Mf = Ta, which has the same structure as the infla- 
tion rule of equation (J5|. Here we assume that either the 
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Fig. 3. Phase difference between two waves originating 
from traversing a layer with index of refraction n 7 and 
thickness d 1 . 



surrounding medium has the same index of refraction as 
medium A or one has to add the corresponding transfer 
matrices Tqa and Tac at the edges, respectively. 

In particular, we are interested in the calculation of the 
transmission coefficient T (also known as transmittance) 
of the light through the stack, which is defined as T = 
\E^\ 2 /\E®\ 2 and can be derived from equation ( p~5] ) by 
eliminating E® . Assuming that there is no incident light 
from the right, i.e., E(* = 0, one obtains 



T = 



Ml 2 + 2detM ' 



(18) 



where \M\ 2 corresponds to the sum of the squares of all 
four matrix elements of M. Further, for the complete stack 
it can be shown that detM = 1, which simplifies the ex- 
pression for the transmission coefficient T. 



3 Results 

In this section we comprise the results for several metallic- 
mean quasicrystals. Especially, we focus on the change of 
the transmission coefficient T in dependence on the used 
inflation rule V and the incidence angle 

For the choice of the wave length of the incident light 
Ao, we use the commonly applied quarter wave length con- 
dition riAdA — rtB&B — Ao/4, which leads for an incidence 
angle d = to an identical optical wave path for the 
light in the two materials A and B of the stack [23 24)25] . 
Hence, equation (14) results in 



A ° Q U Q 
<Z? 7 — — 7TCOS17- = 7TCOS17 

A ujo 



7 ' 



(19) 



In Figure [4] the results for the transmission coefficients 
T(uo) are shown for different approximants of the octonacci 
chain C^ s for an incidence angle $ = 0. Due to the struc- 
ture of the interface and propagation matrices only the 
ratio of the indices of refraction u = ua/tlb is important. 
Here we choose a fixed ratio u = 2/3, and nc = tla for 
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Fig. 4. Transmittance T(uo) for different approximants of the octonacci chain with incidence angle $ = and indices 
of refraction via = tic and u = tia/tlb = 2/3: (a) C^ g with / 5 = 41 layers, (b) C^ g with /7 = 239 layers, (c) C$ g with 
/g = 1393 layers, and (d) C^f with fn =8119 layers. 



the surrounding medium. Results are only displayed in 
the interval [0,cj/cjo] because the pattern of the transmit- 
tance T repeats every uj/ujq for d = (cp. Figure |5j. In 
general, the transmission spectrum of metallic-mean se- 
quences or Cantor sequences is assumed to show either 
complete transmission or complete reflectance for infinite 
systems [19129] similar to the singular continuous energy 
spectrum of quasiperiodic systems [4 6 7]. This character- 
istic becomes more and more visible for larger system sizes 
f m . However, especially for small approximants there still 
occurs a significant transmittance T(uo) > in the pho- 
tonic band gaps since the waves can tunnel through the 
stack to a certain extent. Likewise, T(uS) = 1 is often not 
reached within the bands. 

When varying the incidence angle, we have to keep in 
mind that total reflection occurs for n# < via at the in- 
cidence angles d > $totai = arcsinu -1 and hence we can 
only expect to obtain transmission coefficients for inci- 
dence angles d < #totai- 

The transmittance T($, a;) in dependence on the re- 
duced frequency uj/ujo and the incidence angle $ < $ total 
is shown in Figure [5] and Figure [6] for different metallic- 
mean quasicrystals. Figure [5] displays the transmission for 
the 8th approximant of the octonacci sequence for a broad 
range of frequencies for two examples corresponding to the 
cases ua < tlb and n# > via- As mentioned, for d = 
the transmission spectrum is periodic with respect to the 
frequency. However, for all angles d > the photonic 
transmission bands bend towards higher frequencies for in- 
creasing incidence angles and this effect becomes stronger 
for higher frequencies. Further, for u > 1 one obtains a 
relatively large range of frequencies with nearly complete 
transmission for large incidence angles and for u < 1 there 
are certain frequencies with nearly complete transmission 
over a large range of angles due to the bending of the 
transmission bands. 



In Figure [6] the transmittance T(#, a;) is displayed for 
different metallic-mean sequences. Comparing with the re- 
sults obtained for a periodic alignment of the layers A and 
B in Figure [7J one can see that there is a typical photonic 
band gap at uo = uoq/2 for the periodic chain for d = 0. 
Also for uj/ujo = 6+1/2 (b G No) such gaps occur and 
with increasing incidence angles these photonic band gaps 
bend also towards higher frequencies. We denote the corre- 
sponding frequencies of these photonic band gaps as mid- 
gap frequencies. For the quasiperiodic sequences this band 
gap appears to be much wider, but is intersected by an 
increasing number of new lines of moderately high trans- 
mission for increasing a of the inflation rule. For cu/cuo = b 
(b G No) the characteristics of the transmittance T(#, uo) 
hardly change for the different construction rules. The be- 
havior at the photonic band gap is in consistency with the 
observation that for uj = ujo/2 ($ = 0) the quasiperiod- 
icity is most effective [30 j. For u > 1 the narrow trans- 
mission bands caused by the quasiperiodicity do not bend 
as strongly as the transmission bands already present for 
the periodic case. Therefore, for large incidence angles $ 
the bands corresponding to the periodic case intersect the 
photonic band gaps and complete transmission occurs. In 
this case the transmission spectrum for the quasiperiodic 
and the periodic case look quite similar. For u < 1 all 
bands bend nearly equally strong but also form regions of 
almost complete transmission at certain frequencies for a 
large range of incidence angles (cp. Figure[5|b)). Addition- 
ally, in both cases we obtain that the additional transmis- 
sion bands are more uniformly distributed with increasing 
a of the inflation rule and hence the photonic band gaps 
become smaller. 

The dependency of the transmission bands on the inci- 
dence angle can be explained by the difference of the opti- 
cal light paths in the layers A and B for d > 0. In this case 
we expect that the transmission coefficients T(cj,$ = 0) 
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Fig. 5. Transmittance T(i?, a;) of light through a stack of the materials A and B, which are arranged according to the 
octonacci sequence C$ g with f$ = 577 layers. Results are shown in dependence on the reduced frequency u/uo versus 
the incidence angle d for ua = Tic and (a) u = 3/2 and (b) u = 2/3. 
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Fig. 6. Transmittance T(#, lj) for different quasiperiodic inflation rules: (a) the Fibonacci sequence C^ u with / 14 = 610 
layers, (b) the octonacci sequence C$ g with /g = 577 layers and (c) the bronze mean quasiperiodic sequence Cq y with 
fs = 469 layers. The red lines show the expected limit behavior of the bending according to equation (25). The ratio 
of the indices of refraction is u — 2. 



and T(cj / , $ > 0) yield almost the same results if the over- 
all phase difference for the stack is equal. At first, we con- 
sider only the behavior arising from the change of the 
phase difference within one layer. Hence, for layer B the 
relation 



(p B (u,# = 0) = (PB(w',tiB >0) 



(20) 



has to be fulfilled, which yields a behavior according to 



a/(tf) 



|T9) = 0) (12) 



= 0) 



COS^F 



cos (arcsin (u sin #)) 



(21) 



Note that for tla = tlq the relation $a = $ holds. Like- 
wise, the relation for layer A can be derived with 



|T9} u(d = 0) 



cos$ 



(22) 



However, the actual bending of the transmission bands is 
given as a superposition of equations ([2l|) and ([22]) . For a 



periodic stack with the same number of layers A and B 
the total phase difference ip = (if a + ^ J e)/m/2 results in 
a dependency 



2u;(tf = 0) 



cos d + cos (arcsin sin $)) 



(23) 
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Fig. 7. Transmittance T($, cj) for a periodic sequence of alternating layers ABAB . . . with length 500 for the ratios 
of the indices of refraction (a) u = 2 as in Figure [6] and (b) u = 0.5. The additional lines in (a) and (b) show the 



expected limit behavior of the bending: the green, blue, and red lines correspond to the equations (21), (22), and (23) 



respectively. Plot (c) shows these limit functions as well but additionally the bending for quasiperiodic sequences in 
dependency on the metallic mean r according to equation (25). The Moire patterns in (a) and (b) are caused by strong 



oscillations of the transmittance T, which are not adequately displayed due to the used resolution. 



The functions (|2l|)-(23) are shown in comparison with the 
numerical results of the transmission in a periodic stack 
of layers A and B in Figure a) and (b). The functions 
(23) fit the bending of the transmission bands quite well, 



especially for small angels. Further, we obtained that the 
regions of almost total transmission behave according to 
equation ( [21] ), which indicates that cps plays a significant 
part for this behavior. In contradistinction the bending 
caused by cpA often coincides with regions of lower trans- 
mission coefficients. This effect becomes even more evident 
for systems with w <C 1 or u > 1. 

For the quasiperiodic stacks the numbers of layers A 
and B are not identical but related to the metallic means 
r by equation Q . The overall phase difference in the limit 
of infinite stacks is then given by 



1 



1 



-<Pb fi 



and yields a bending according to 



(1 + t)u;(i? = 0) 



3 $ + cos (arcsin (u sin #)) 



(24) 



(25) 



These functions are shown in Figure [7^c) for the two cases 
u = 0.5 and u — 0.2. Again, we obtain a bending in be- 
tween the behavior arising from one single layer of medium 
A or B. In particular, the bending increases with the pa- 
rameter r (respectively a) in the case u < 1 and shows 
the reverse behavior for u > 1. Further, for the quasiperi- 
odic systems the bending according to equation (25) ap- 
proaches for r — » 1 the bending of a periodic stack and 
for r ^> 1 the behavior arising from one single layer B. 
Comparing the functions (25) with the numerical results 



of the transmission for the quasiperiodic layered systems 
with different values of r (cp. Figure [6]), we find a very 
good resemblance of the bending behavior for small and 
intermediate angles. 

In Figure [8] the characteristics of the transmittance T 
are investigated with respect to the change of the ratio 
of the indices of refraction u. In Figure a) -(c) we study 
the influence of the reduced frequency uj/ujq and in Figure 
[8jd)-(f) the change of the transmittance T with respect 
to the incidence angle $ for the periodic, the golden mean 
and the silver mean sequence. 

In the first case, i.e., for the variation of the frequency, 
it becomes again clearly visible that the quasiperiodicity 
has the largest effect for uj w cjo/2. In this region the 
transmittance T(uo,u) strongly varies with the used con- 
struction rule. By increasing the parameter a, more and 
more bands appear. In the periodic case complete trans- 
mission occurs for this frequency only for u = 1 , whereas 
for the quasiperiodic systems there are several bands in 
this region with a high transmission coefficient. Of course, 
complete transmission occurs always for u = 1. 

The results for the investigation of the influence of the 
incidence angle $ for different refraction coefficients are 
shown in Figure |8fd) -(f). A comparison of the plots for 
the periodic and the metallic-mean sequences shows that 
there is a similar structure of the transmittance T($, u) 
near the regions of total reflection. However, for small in- 
cidence angles there is a completely different behavior, i.e., 
we have zero transmission for u ^ 1 in the periodic case, 
whereas for the quasiperiodic sequences complete trans- 
mission occurs for a significant range of ratios u. 
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Fig. 8. Comparison of the transmittance T of light through a stack made of the materials A and B for the periodic, 
golden, and silver mean systems, i.e., the same systems as in Figure [7^, |6^i and ^jp. Results are shown in panels (a)- 
(c) in dependence on the reduced frequency uj/ujq and the ratio of the refraction indices u = ua/tlb for a constant 
incidence angle d = 0.2 and in panels (cl)-(f) in dependency on u and d for uo = ujq/2. The white regions in the lower 
right corners in (d)-(f) correspond to cases of total reflection. 



4 Conclusion 

We investigated the transmission of light through a stack 
of quasiperiodically aligned materials with respect to the 
underlying inflation rule and for different incidence angles. 
One main result for the quasiperiodic systems is that mod- 
erately high or even complete transmission can occur for 
certain ranges of parameters, for which no transmission is 
found for the periodic system: e.g. we obtained additional 
bands of moderate transmission for frequencies near the 
midgap frequency (b + 1/2)ujq (b G No) with the number 
of bands increasing depending on the parameter a of the 
inflation rule, and for a whole range of ratios of refraction 
indices u even complete transmission occurs. 

Further, for all considered systems almost complete 
transmission occurs for incidence angles d close to the an- 
gle of total reflection. This is caused by the bending of the 
transmission bands due to the different optical light paths 
in mediums A and B. Additionally, the quasiperiodicity 
shows only an influence on the transmission coefficients 
around the region of the midgap frequencies (b + 1/2)cj 05 
whereas in between the patterns appear to be similar for 
different quasiperiodic generation rules as well as for the 
periodic system. 
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